Microtubule Dynamics and Oscillating State for Mitotic Spindle 
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We present a physical mechanism that can cause the mitotic spindle to oscillate. The driving 
force for this mechanism emerges from the polymerization of astral microtubules interacting with 
the cell cortex. We show that Brownian ratchet model for growing microtubules reaching the cell 
cortex, mediate an effective mass to the spindle body and therefore force it to oscillate. We compare 
the predictions of this mechanism with the previous mechanisms which were based on the effects of 
motor proteins. Finally we combine the effects of microtubules polymerization and motor proteins, 
and present the detailed phase diagram for possible oscillating states. 
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I. INTRODUCTION 

Establishment of geometrical polarity during mitosis 
which is a result of asymmetric cell division is a funda- 
mental fact in many live systems. The understanding 
of underlying molecular mechanisms for this phenomena 
is an important and challenging issue in biophysics |l|. 
Including cell division, various functionalities of cell are 
related to the microtubules, the most rigid filaments of 
eukaryotic cytoskeleton. These filaments are responsible 
for the mechanical properties of cells @, y[ • The func- 
tional properties of microtubules are determined by the 
structures of their monomers and the interaction mecha- 
nism packing them into rigid polymers. These long and 
hollow cylinders are made from a — /3 tubulins, proteins 
in the eukaryotic cells [4| . 

Microtubules growing from two organizing centers in- 
side the cell, form a reliable scaffold for the mitotic spin- 
dle and make the mechanics of the cell. In the ani- 
mal cells during mitosis, microtubules are nucleated from 
centrosomes, organizing centers near nucleus [2|. Some 
microtubules grow from one centrosome toward another 
and construct a firm bridge between two centrosomes. 
These microtubules with two centrosomes, form the com- 
plicated structure of the mitotic spindle which plays an 
important role in cell division. Microtubules also, are 
required for divorcing the chromosomes at the end of mi- 
tosis [5[. Some other microtubules, astral microtubules, 
growing from centrosomes, reach the cell cortex and push 
the spindle body. They are absolutely required for cor- 
rect positioning of mitotic spindle |6| . 

It is understood that the spindle motion during the 
cell division is responsible for asymmetric cell division. 
To gain insight into the cellular mechanisms by which 
the spindle oscillates and repositioncs during the mito- 
sis, extensive studies are carried out on the C. Elegans 
embryo |7| . All experiments have clarified the importance 
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of microtubules and their polymerization. 

Considering a solid spindle body and astral mi- 
crotubules, interacting with cell cortex, two different 
deriving mechanisms for spindle motion are possible. 
The first is the pushing and pulling forces related to the 
motor proteins connecting the microtubule tips to the 
cell cortex and the second is the pushing mechanism due 
to the microtubule polymerization. The first mechanism 
is studied by S. W. Grill et al. |8J], where, they have 
shown that the resulting forces make the spindle oscil- 
lation possible. Regarding the second mechanism, the 
existence of pushing forces is experimentally investigated 
and theoretically studied [9l-lll|. In such systems the 
effects of microtubule bending and buckling arc discussed 

In this article we concentrate on spindle body motion 
with deriving forces originated from the microtubule and 
cell cortex interaction. We present a simplified theoreti- 
cal model with minimum requirement to generate an sta- 
ble oscillating state for spindle body. In section II, we 
briefly review the governing equations for the dynami- 
cal instability of microtubules and also Brownian ratchet 
model[15|,[l6[. In section III, we show that the interaction 
of microtubules with cell cortex can produce an overall 
effective mass for the spindle body and force it to oscil- 
late. The effect of the molecular motors is considered in 
section IV. Discussion about our results and concluding 
remarks are presented in section V. 



II. MICROTUBULE'S DYNAMICS AND 
BROWNIAN RATCHET MECHANISM 

Let us consider a bundle of stiff and noninteracting 
microtubules growing from a centrosome (nucleating cen- 
ter). The life of these microtubules can be seen as two 
different phases: the growth and the shrinkage phases. 
To describe the stochastic transition of microtubules be- 
tween two phases in a simplified model, four parame- 
ters are defined. Two transition rates, f cat and f res , 
the catastrophe and rescue frequencies that quantify the 



rates by which the microtubules change their growing and 
shrinking phases and also average growth and shrinkage 
velocities for microtubules at the bulk which are shown 
by Wg and v s respectively. Denoting by n + (x, £), the num- 
ber density of growing microtubules at position x and by 
n~(x, t), the number density of shrinking microtubules, 
the d yna mical equations can be written in the following 
form (171 ]: 



dt 
d 

df 



n (x, t) 
n~(x, t) 



■Sn + (x, t) - f C atn + {x, t) + f res n (x, t), 

■Sn~(x,t) - fresTl~(x,t) + f C at11 + (x,t), 



where Sn = n (x ^f l,t) — n {x, t) and I stands for the 
monomer length. 

Despite the above equations for microtubule dynam- 
ics in the bulk, the microtubules have different evolving 
equations on the centrosome (nucleation cite), and also 
on the cell cortex. To express the dynamical equations 
for the microtubules on the centrosome at x = 0, we de- 
fine v as nucleation rate . The nucleation rate is a rate 
by which free inactive sites on the centrosome change to 
active sites and grow by absorbing tubulins. Now the 
boundary equations at the centrosome read [181 ]: 



* B+M 



vn-(0,t)--e-n + (p,t), 



j t n~(0,t) = -vn-(p,t) + jn-(l,t), (1) 

Correspondingly the boundary equations at the position 
of the cell cortex that is located at a fixed position x — L, 
can be written as (la: 



where A 



/re. 



The growth and shrinkage velocities that we have used 
in above descriptions, are assumed to be the velocities 
for free microtubules in the bulk. For a microtubule 
tip reaching an obstacle, polymerization process can per- 
sist by pushing which repositions centrosome (nucleation 
cite). Brownian ratchet model is a physical mechanism 
that describes the dynamics of this process [l5|, [la, [l9| . 
Denoting by F, the force that the wall exerts on the tip 
of a microtubule, the following linear response, force ve- 
locity relation can be written: 



F 



Fstall 1 

V a 



(6) 



where v is the polymerization velocity in the presence 
of the load and v g stands for the polymerization veloc- 
ity in the absence of the load. The stall force F sta u, is 
the threshold force that stops the polymerization pro- 

k B T (-, _ k °ff\ 



cess. The stall force is given by F sta ii 



where k® n and h®** are polymerization and depolymer- 
ization rates. One should note that the bulk polymeriza- 
tion velocity is related to the microscopic rates by: 



' (%n 



k off) 



(7) 



Defining friction coefficient by £o 
the force velocity as: F — £o (v g - 



F 3ta 



we can rewrite 



v). 



In the next section we will use above description about 
microtubule dynamics and Brownian ratchet mechanism 
to analyze the dynamics of spindle body. 
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■n + {L -l,t)- f wa iin + (L,t), 
■^-n-(L,t) + f wall n + {L,t), (2) 



where f wa ii represents the catastrophe frequency on the 
obstacle. 

For the next step, we present the steady state solu- 
tions to the above system of equations. Using matrix 
techniques for difference equations and for a wall that is 
fixed at position L = Lq, the solutions read 18 |: 



n + (x) = Aexp(— — ) 
A 



n (x) 



v, 



.r , 



Acxp(--) 



v s A' 

and the boundary values are: 



< x < Lq 
< x < L , 



(3) 



n"(0) 



l'n 



'-n+(0), n+(L ) 



Jwai 



-n-(Lo), (4) 



with 



A= [ A(l + ^)(1 
v» 



-^\ i u g 



V 



Jwali 



(5) 



III. DYNAMICS OF SPINDLE BODY 

Astral microtubules with the mitotic spindle constitute 
our one dimensional model shown in Fig.[T] In our model, 
spindle structure, behaves like a rigid body. Dynamical 
instability of astral microtubules makes stochastic tran- 
sition between growth and shrinkage phases. During the 
polymerization process, each microtubule that reaches 
the wall, pushes it with force F that obey the Brownian 
ratchet response function [l5j . Additionally we assume 
that the cortex is a fixed and rigid wall that does not 
displace or deform during spindle motion. The mitotic 
spindle then, feels the reaction of microtubule forces on 
both sides. The reaction force from the left (right) hand 
side wall is proportional to the number of microtubules 
which have reached the left (right) wall. Denoting the 
number of attached microtubules in the left side by Ni 
and the corresponding number at right side by N ri and 
using the linearized Brownian ratchet response function, 
the net polymerization force acting on the spindle body 
can be written as: 



-t ii 



(Ni-NJtoVg-iNt+NJZoL. 



(8) 
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FIG. 1: (color online). One dimensional modeling of spindle 
body and astral microtubules in mitosis; two centrosomes and 
the microtubules between them assumed to be rigid and form 
spindle. The cell cortex is fixed and the net force due to the 
pushing of attached microtubules cause the spindle to move. 



can be expanded as below: 

Ni = N? + Nfe™*, N r = N° + A^V"*. (12) 

We note that the time scale for the polymerization pro- 
cess is smaller enough than the time scale for the spindle 
body motion. This fact allows us to assume that for a def- 
inite value of L(t), the number of attached microtubules 
is related to the length L(t) through the steady state re- 
lations presented in Eq. 01 We also expand the number 
density up to the first order of perturbation, for zeroth or- 
der value, we can use the time independent steady value 
(i.e. no oscillation). Now, to find number densities on 
cortex, we assume: 



which L is the spindle's velocity. Beside the above force, 
the effects due to the complex properties of cellular fluid 
should be added. The viscoelastic properties of the cy- 
toplasmic fluid can be modeled by defining an effective 
viscous drag coefficient (i.e. £) and an effective clastic 
modulus for the spindle body (i.e. n). The effects of 
other astral microtubules in the real three dimensional 
configuration are also included in this elastic constant, 
n. As a result, the governing dynamical equation for the 
spindle body reads: 



£L + k(L - L Q ) = F, 



pol • 



(9) 



The elastic term acts as a centering force, it enforces the 
spindle body to choose an equilibrium position Lq. Due 
to the dynamic instability of microtubules structure, a 
growing microtubule may change to a shrinking one even 
if it has reached the cortex surface. However, this hap- 
pens with a new catastrophe frequency, compared to its 
bulk value. Then, the number of attached microtubules 
at both sides obey the following equations: 



Ni = -UallNi + (v g - L)Nn+(L,t), 
N r = -fwallNr + (v g + L)Nn+(L,t). 



(10) 



which N is the total number of microtubules. In the 
limiting case, where the spindle body doesn't move, these 
equations are the same as Eq. [2] multiplied by N, as a 
total number of microtubules . 

To solve above dynamical equations and consequently 
obtain the dynamical behavior of spindle body, we as- 
sume that the system eventually reaches an oscillating 
steady state that is a steady state solution which oscil- 
lates around an equilibrium state. We denote the equi- 
librium value of spindle position by Lq and the number 
of attached polymers by N® and N®. Assuming that 
the polymerization forces in comparison with the clastic 
restoring forces are small, we are looking for perturbativc 
solutions for spindle body in the following form: 



L = L Q + L ie " 



(11) 



where L\ is a small deviation from the equilibrium posi- 
tion. In this case the number of attached microtubules 



n+(L,t)^n+(L )+L 1 e luJt —n+(x) x=Lo , (13) 



nt(L,t)=nt{L )-L 1 e iut —rf(x) x = Lo , (14) 
then we will have: 



d 
dx 



n+(L,t) - n+(L,t) = -2^n+ (L )e iut , (15) 



and 



~(L,t) + n+(L,t) = 2n+(L ). 



(16) 



If we put the above information in Eq.® and Eq. ([T0|) . 
we may derive F po i which can be used in Eq.(j9|). In the 
first order of perturbation expansion, our governing dy- 
namical equation will be: 



(Xo - *XiW + X2U 2 ) L x = 0. 



(17) 



As one can see the overall motion of the spindle body is 
determined by a viscose term xi and an inertial term \2- 
These effective parameters are given by: 

X0 = f W all(v g \B - j), xi = fwallO- + 2B) + |(18) 



X2 = 1 + B, B = 



2Z a v g Nn + (L ) 



SsJwall 



(19) 



To investigate the existence of any oscillating state with 
nonzero amplitude, we can solve the above equation and 
find the allowed values for lu. Separating the real and 
imaginary parts of the frequency as lu = u r + iwj, we can 
study the motion for different values of parameters. 

Fig. [21 and Fig. [3] show two different phase diagrams 
of the system. Parameter values which we have used 
are: v g = lAl/im/min, v s = 30.7 ^irn/min, f res = 
0.015s" 1 , f cat = 0.0037s" 1 , k = 4 x 10" 6 A^/m[|, 
£ = 10~ 6 N s/m®, KbT = i.lpN.nm, L = 10/xm, 
v = 0.01s" 1 , I = 0.6nm and k on = 100s" 1 H3|. The 
phase diagram according to parameters f wa ii and f ca t 
(Fig. [2]), shows that for catastrophe frequency on the wall 
(fwa,u), below a critical value, spindle can oscillate. Fig. [3] 
shows the phase diagram for different values of the num- 
ber of microtubules (N) and catastrophe frequency on 
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FIG. 4: (color online). One dimensional modeling of the spin- 
dle, microtubules and motor proteins. Motor proteins on the 
cell cortex can attach to the microtubules and exert force on 
spindle body. Microtubules polymerization also exerts force 
on the spindle and these two effects cause spindle to move. 



FIG. 2: (color online). The phase diagram for different values 
of fwaii and feat for a given experiment [20| . Three different 
states are separated. For states with uj r — 0, the oscillating 
motion is not possible and all solutions are damped to zero, 
while for oj r =^ 0, oscillating solutions are possible. For ui r > 
uii, the damping time scale is large enough to have a well 
defined oscillation. 
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FIG. 3: (color online). The phase diagram for different val- 
ues of total number of microtubules N and f wa ii for a given 
experiment [2y|. Three different states are separated. For 
states with u: r — 0, the oscillating motion is not possible and 
all solutions are damped to zero, while for uj r =^ 0, oscillating 
solutions are possible. For oj r > uii, the damping time scale 
is large enough to have a well defined oscillation. 



decrease to 0.0018 s _1 . Clearly by changing the unknown 
parameters, we can tune the frequency to capture the ex- 
act oscillation frequency observed in live systems [8j. 



IV. MICROTUBULES, MOTOR PROTEINS 
AND SPINDLE 



In this part we will consider the effect of motor pro- 
teins [8| as well as microtubules on the spindle's motion. 
As shown in Fig. [4j microtubules grow from the centro- 
somc toward the wall and when they reach the wall, they 
exert force on it. Motor proteins on the cortex can bind 
to these microtubules and exert force on them as well. 
We assume that some fraction of motor proteins on the 
cortex, can attach to those microtubules which are not 
farther away than a critical distance from the cortex. In 
this model we consider both forces that are exerted from 
microtubules and motor proteins and try to find the over- 
all equation of motion for the spindle body. 

We first present a simplified and mean field version of 
the motor protein mechanism introduced by S. Grill et 
al. [8j. A motor protein can attach to a filament, walks 
along it and exerts force on it. Denoting the walking 
velocity by v and corresponding force by / and for small 
velocities, the linear response: v = vq (1 — f/fo) holds 
during this motion [8|, [l3| . Here /o is the stall force by 
which the motor protein stops and vo is the velocity of 
motor in the absence of any force. Denoting the elastic 
properties of motor protein with a single parameter k and 
its extension by x, we can write the force motor feels as 
/ = kx (FigHJ). Then the motors velocities on the right 
wall and left one become: 



the wall (f W aii) ■ It can be seen from this graph that 
there is a critical value for N (depending on the value of 
fwaii), below it the oscillatory state is not present. 

As an example for the numerical value of the spindle 
oscillation frequency (Fig. [21 Fig. [3]), and for numer- 
ical values: N = 100, f waU = 0.001s" 1 and f cat = 
0.007s _1 , the oscillating frequency for the systems is: 
u> r = 0.0023 s _1 . Changing the wall catastrophe rate 
to feat = 0.005 s" 1 , we sec that the oscillation frequency 



xl = v - hxl + L, 
xr = v - fix R - L, 



(20) 
(21) 



where L is the spindle velocity, xl is the extension of 
the linker (spring which connects the motor protein to 
the cortex) on the left side, xr is the extension of the 
linker on the right side and fi is equal to kvo/fo. The 
motor protein can attach to microtubule with a rate oj on 
and detach from it with a rate ui ff- According to these 
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FIG. 5: (color online). Phase diagram of possible states for 
different values of the rate u) on (in unit of s" 1 ), and motor pro- 
teins number M . The colored region bounded by solid lines, 
shows the states of stable oscillatory solutions for the case 
where the effects of motor proteins and polymerization forces 
are included. Dashed line shows the boundary for oscillatory 
and non oscillatory states when the polymerization force is 
not included yj- Parameter values used in this graph are: 
k = 8.3 x l(T 3 A7m, f = 3pN, v = 1.8/im/s , w = 5s" 1 0, 



L = 10/xm, N = 100, f wa ii = 0.001s" 



= 1. 41 fim/min, 



30.7 fim/min, f res = 0.015s -1 and f cat = 0.0037s" 1 [17 



rates, we can find the fraction of bond motor proteins as: 

QbL = -^offQbL + Won(l ~ QbL,), (22) 

QbR = —UoffQbR + Won(l - Q&fl)) (23) 

QbL is the fraction of bound motors on the left side while 
QbR is the fraction of bound motors on the right side. 
For simplicity we have ignored the diffusion of motor 
proteins on the filament. We can replace all time scales 
corresponding to the walking processes on the filament 
by w~ft, the average time that a motor spends in at- 
tached state. Since the detachment rate (i.e. w //) de- 
pends on the force motor feels, it depends on the linker's 
length. However the attachment rate (i.e. uj on ) is con- 
trolled by the barrier that motor feels to connect to the 
filament and local temperature, and hardly depends on 
linker's length^]. The force dependent detachment rate 
is co ff = Wo cxp(ka\x\/ ksT) where wo is the detachment 
rate for motor proteins in the absence of any force, k is 
the spring's constant, x is the linker's length and a is the 
molecular length scale. @ 

Two crucial elements in Grill et al. mechanism are 
the stochastic nature of polymerization process and the 
force dependent detachment rate that result oscillatory 
motion for the spindle. Keeping these two elements, we 
also have considered the force due to the microtubules 
polymerization. If microtubules reach a critical distance 
from the cortex, motor proteins have the chance to attach 
to them with the rate w on . If we denote by Pl and P R , 
the number of such microtubules on the left and right 
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FIG. 6: (color online). In this phase diagram we compare 
the results of a model that is based on the polymerization 
forces and a model that takes into account both the poly- 
merization forces and the effects due to the motor proteins. 
As in Fig. [3] the dotted line shows the boundary of differ- 
ent phases for the case where only the polymerization effects 
are considered. Including the effects due to the motor pro- 
teins, we see that the dotted line changes to dashes-dotted- 
dotted line and also a new phase of stable oscillating state 
(i.e. uji < 0) comes to existence. This new sate is shown as 
a colored region. For the physical parameters we have used 
the following values: v g — 1.41fj,m/min, v s = 30. 7/im/mm, 



fres = 0.015S- 1 , fcat = 0.00373" 

/o = 3pN, v - 1.8/oti/s , ujo 
Lo = 10^m and M = 100. 



n/rm 

m 

5s" 



k = 



8.3 x 10~ 6 N/m, 
n = 0.25s" 1 B, 



sides, we can write: 

p L= (nt(x,t) + n£(x,t))dx, 

Jl-i 

Pr= (n^ i (x,t) + n] i (x,t))dx, 

Jl-i 



(24) 

(25) 
(26) 



where Iq is the critical distance from the cortex for a mi- 
crotubule by which motor proteins can attach to it. Here 
nj j {x,t) (rift(x,t)) is the number density of growing mi- 
crotubules on the left (right) side and n^(x, t)(n^(x, t)) 
is also the number density of shrinking microtubules on 
the left (right) side. 

Combining the effects due to the polymerization and 
also the motor proteins, we can arrive at the following 
dynamical equation for the spindle body: 

£Z + k(L - L ) = (N L - N R )£ v g - (N L + N R )£ L 

+M R P R Q bR kx R - M L P L Q bL kx L ,(27) 

where Nl and N R are the total number of microtubules 
in the left and right sides and Ml and M R are the total 
number of motor proteins in the left and right sides. 

Now we can solve Eq.fTUland Eq. I201I271 and obtain the 
overall dynamical properties of the spindle body. Follow- 



ing the same perturbative method described in the pre- 
vious section, we can arrive at the following equation in 
Fourier space: 

X (w) xii=0, (28) 

where the dispersion x(w) is given by: 

XM = X4^ 4 + X3^ 3 + X2W 2 + xiW + Xo, (29) 

and the coefficients are given by: 



X4 

X3 

X2 
Xl 

s 

D 
F 

H 



1 + B, xo = »C{Ef wall - v g XB), 
(l + B)(H + f waH ) + B + E + F, 
(1 + B)(jtC + UaiiH) + E(H + f wall ) 
-v g \B + BH + F(f waM +C)-G, 
(1 + B)fiCf wall + E(fiC + Hf wall ) 
-v g \BH + liCB + FBf wall - Gf wall , 
2^ v g Nn+(L ) 



t,Jwai 



C = Won +V ff(Xo), 



^, E = 2MD^XP 0+ « 

2kMP C kv , 
-, , G = w {x )F, 



\x + C, fx 



kvp 



To analyze the possible states we investigate the numer- 
ical solutions to x( w r + *£*>$) = 0. Fig.[5j shows the phase 
diagram of different states of the spindle for different val- 
ues of M, number of motor proteins, and w on , attachment 
rate for motor proteins. There is a range of values for ui on 
and M that spindle has stable oscillations. This graph 
corroborates our main idea that microtubules polymer- 
ization is not negligible and changes the quantitative be- 
havior of oscillation. 

In Fig. |H1 we have shown the phase diagram of the sys- 
tem for different values of f wa ii, catastrophe frequency 
on the wall, and N, number of microtubules. In this 
graph, we compare the results of a model that is based 
on the polymerization forces and a model that takes into 
account both the polymerization forces and the effects 
due to the motor proteins. As in Fig. [31 the dotted line 



shows the boundary of different phases for the case where 
only the polymerization effects are present. Including 
the effects due to the motor proteins, we see that the 
dotted line changes to the dashed-dottcd-dotted line and 
a new phase of stable oscillating states(colored region) 
comes to existence. As an example of the oscillating 
frequency for f waU = 0.0006 s~\ M = 100, N = 100 
and ui on = 0.25 s" 1 we find that u> r = 0.006 s^ 1 and 
Ui = -0.0005 s" 1 . 



V. CONCLUSION 

Similar to the Grill et aVs mechanism, we believe 
that the stochastic nature of the polymerization process 
and also the force dependent polymerization velocity are 
two important points that derive the oscillations in our 
model. The effects due to the motor proteins are also 
considered in the frame work of a simplified version of 
the mechanism presented in Q. We have added both 
the polymerization forces and motor forces to obtain a 
more realistic description of the spindle motion. Taking 
into account these two mechanisms, we have studied the 
phase diagram for spindle dynamics. We have presented 
detailed phase diagrams for spindle motion and compared 
the phase diagram with and without the microtubules 
effect (Figf5]). We see that adding microtubules poly- 
merization effect changes the phase diagram; the critical 
number of motor proteins for which oscillation starts, are 
different in two cases. There is also a critical number for 
microtubules, above which the oscillation occurs. This 
critical amount changes if we consider/ignore effect of 
motor proteins. The numerical values for the oscillation 
frequency are in the range of observed experimental val- 
ues. 

In conclusion, while our results confirm Grill ei.a/.[8( 
former idea that motor proteins are the main source of 
spindle oscillation, we see that microtubules polymeriza- 
tion has non-negligible effect on the quantitative behavior 
of such oscillation and can not be ignored. 
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